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1. Introduction 

Non-zero quadrupole strength in the N to A transition is interpreted as a signal for deformation 
in the nucleon and/or A. Hadron deformation is a fundamental property that depends on QCD 
dynamics and determining the shape of the nucleon has motivated experimental measurements of 
quadrupole form factors in the N to A transition for many years. Selection rules in the N to A 
transition allow for a magnetic dipole, Ml, an electric quadrupole, E2, and a Coulomb quadrupole, 
C2, with the corresponding transition form factors Gm\, Ge2 and Gqi as functions of the four- 
momentum transfer squared, q 2 . The magnetic transition, which proceeds through a quark spin- 
flip, is expected to be the dominant one while non-zero quadrupole strengths can be obtained in 
many models by assuming a deformed nucleon and/or A as for example in the quark model where 
a D-wave admixture is assumed in the nucleon and A wave functions. Pion electro- and photo- 
production experiments near the A(1232) resonance have clearly established non-zero quadrupole 
strength in the range of a few percent with respect to the dominant Ml strength providing support 
for deformation In ref. [||], a precise measurement at — q 2 = 0.127GeV 2 extracted the values 

Rem = -^- = (-2.3 ±0.3 ±0.6)%, R SM = - P^^ 1 - = (-6.1 ±0.2 ±0.5) % (1.1) 
Umi Zm^ Umi 

for the ratio of the electric and coulomb quadrupole to the magnetic dipole form factors, Rem (or 
EMR) and R$m (or CMR). The first uncertainty given in Eq. ( |ET| ) is the statistical and systematic 
error and the second the model error. At the same time experiments at Jefferson Lab have studied 
the q 2 dependence of these three form factors in the range of 0.4 — 1.8 GeV 2 

Within lattice QCD one can evaluate these transition form factors starting directly from the 
QCD Lagrangian. The first lattice study of the N to A transition was carried out very early 
on [Q] and, although the quadrupole form factors were not determined to high enough accuracy, 
the methodology for the calculation of this matrix element on the lattice was presented. Using 
the formalism laid out in Ref. [Q] a more detailed study followed in Ref. [^] which established 
the existence of small and negative electric quadrupole strength for a given momentum transfer q 
within quenched QCD. In subsequent work ^] we focus on the determination of the momentum 
dependence of the transition form factors. We carry out a number of improvements that enables us 
to extract both electric and coulomb quadrupole form factors to high enough accuracy to exclude a 
zero value. The improvements come from optimizing the interpolating field used for the A so that 
the three-point function is evaluated for the maximum allowed values of lattice momentum vectors 
for a given q 2 and from using a simultaneous overconstrained analysis of the resulting matrix ele- 
ments. We obtain small negative values for the ratios Rem and R$m in qualitative agreement with 
the experimental measurements for values of — q 2 in the range 0.4 to 1.5 GeV 2 . On the other hand, 
the values obtained for R$m are smaller than experiment for — q 2 < 0.4 GeV 2 . Similarly the values 
obtained for Gmi in this quenched study, after linearly extrapolated to the chiral limit, disagree with 
the experimental measurements raising questions about the importance of sea quark effects and the 
validity of the linear extrapolation of the lattice data. Recently a relativistic chiral effective-field 
theory calculation of the Rem and R$m ratios shows that the linear extrapolation is incorrect at 
low q 2 and a stronger dependence on the pion mass sets in at a pion mass value where the A be- 
comes unstable. This stronger mass dependence drives R$m to more negative values resolving the 
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discrepancy between lattice results and experiment at low q 2 . This study explicitly demonstrates 
the importance of pion cloud contributions in the determination of these transition form factors. 

In this work we attempt to study the effects of unquenching on the three transition form fac- 
tors using dynamical MILC configurations for the sea quarks and employing domain wall valence 
fermions. Given the fact that chiral fermions are still very expensive to simulate this "hybrid" 
scheme appears as a reasonable compromise and has already been employed in heavy quark spec- 
troscopy and generalized parton distributions calculations 0. 



2. Lattice techniques 



The matrix element for the N to A electromagnetic transition for on-shell nucleon and A states 
and real or virtual photons has the form JT0| ] 



<A(pV) \U\N(p,s)) 



niA mjv 



3 Ga(p') E N (p] 



1/2 



u z (p' ,s')ff T ^u(p,s) 



(2.1) 



where p,s and p',s' denote initial and final momenta and spins and u T (p' ,s') is a spin-vector in 
the Rarita-Schwinger formalism. The operator €? TM can be decomposed in terms of the Sachs form 
factors as 

= G m (q 2 )KZ + G E2 (q 2 )Kl» + G c2 {q 2 )K% , (2.2) 

where the exact expressions for the kinematical functions K xti can be found in ref. [0]. The extrac- 
tion of the Sachs form factors requires the computation of the three-point function G^ N (t2,h;p ',p;T) 
along with the nucleon and A two-point functions G NN and Gjj. The nucleon source is taken at 
time zero, the photon is absorbed by a quark at a later time t\ and the A sink is at an even later time 
?2- Provided the Euclidean time separations t\ and t% — 1\ are large enough, the time dependence 
and field renormalization constants will cancel in the ratio 



R a (t 2 ,ti;p',p;r;n) = 



(Gf( f2 ,p';r 4 )) 



(G NN (t 2 -t u p;T 4 )) (Gf (t ljP ';r 4 )> (Gf (fe,p';r 4 )> 



1/2 



k-»l»l,»l»l 



n ff (p',p;r;M), (2.3) 



(Gf (fe-* 1)P ';r 4 )> <G^( fl ,p;r4)) (G™ (t 2 ,p;U)) 
where a denotes the vector index of the A field. The matrices T are projections for the Dirac indices 

(2.4) 

L \ U U / 

We fix the A at rest and therefore q = p' — p = —p. Q 2 = —q 2 is the Euclidean momentum trans- 
fer squared. Determining n CT (q ;r;/i) for given values of a and T by fitting to the plateau of 
Ro(t2,ti',p\p ;H ;U) enables us to obtain the Sachs form factors. For example, 





n ff (q ;r 4 ;At) = iAe a W p J ' G M x{Q 



(2.5) 



with A a kinematical coefficient. Other similar relations can be found in Ref. [Q]. A novel feature 
of Ref. [^] is the choice of linear combination of ratios R a {t% , t\ ; p p ; F; y.) for appropriate values 
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of (<7,r) such that, for given Q 2 , a maximum number of photon momentum transfers q give non- 
vanishing contributions. The optimal combination for Gmi, referred to as Si -type, is 

3 

Si (q; M) = £ n ff (q ;r 4 ; At) = iA [(p 2 - p 3 )S hil + (p 3 - p^S^ + (pi - Pt)^] G m (Q 2 ) (2.6) 

(7=1 

where all the vectors q related through lattice rotations participate. Similarly, the quadrupole form 
factors are obtained from type S 2 given by 

3 

S 2 (q;n) = £ Uo(q;T k ;n) 

= flj [(p2+P3)5l )M + (p3 +Pl)&2,ii + (pi+Pl)^) G E 2(Q 2 ) 
~ 2^-(p l p 2 +PlP3+P2P3) 

P 

for /l = 1,2, 3. For /i = 4 we have 

6 C 

S 2 (q;n=4) = ^( Pl p2 + pm+P2P3)Gc2{Q 2 ), (2.8) 

where B and C are kinematical factors. The three-point functions needed for the evaluation of Si 
and S 2 are obtained by fixing the appropriate combination of interpolating fields for the A and re- 
quire only one sequential inversion for each fixed sink type. Therefore, two independent sequential 
inversions suffice for the complete evaluation of the form factors for all available lattice momenta 
q. The full set of lattice measurements for any available pL and q for a given Q 2 are simultaneously 
analyzed and the form factors are extracted from a global ^-minimization procedure [Q]. 

The dynamical configurations that we are using are generated by the MILC collaboration and 
are available from the NERSC archive. For the lattice of volume 20 3 x 64 we use the ensembles 
with the strange quark mass fixed at am s = 0.05 and the light quark flavors at am u c i = 0.03 and 
0.02 simulated with the improved staggered (Asqtad) action which ensures better scaling properties 
at the lattice spacing of a = 0.125 fm used in the simulations. In addition we use configurations 
generated on a lattice of size 28 3 x 64 i.e. physical spatial volume of (3.5 fm) 3 at am s = 0.05 and 
amu,d = 0.01. We apply hypercubic (HYP) smearing to the configurations and use Domain wall 
fermions in the valence sector. The size of the fifth dimension is set to L 5 = 16 where it has been 
shown [^] that, for these lattices, the residual quark mass, computed from the divergence of the 
five-dimensional axial current, is at least an order of magnitude smaller than the bare mass. The 
height of the domain wall parameter is set at am$ = 1.7 and Dirichlet boundary conditions are 
imposed at t/a = 32 since only half the time extent of the lattice is needed. Smeared quark wave 
functions are known to increase the overlap of the interpolating fields with the hadronic states 
reducing the contamination from excited states. Before the calculation of the valence propagators 
we apply gauge invariant Wuppertal smearing to the local quark fields q(x) via g sme ar(x, = Ly (1 + 
aH)"(x,y)q(y,t) with H the hopping matrix, a ~ 3 and n = 30. 

The bare Domain wall quark mass parameter (am q ) DW has been tuned by the LHP collabora- 
tion [^] by requiring that the pion mass computed with the Domain wall valence quarks equals the 



G E 2(Q 2 ) + ^^G C2 (Q 2 ) 



(2.7) 
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mass of the lightest pion computed with the staggered fermions. The resulting quark mass param- 
eters and corresponding pion masses obtained by the LHP collaboration are shown in Table 1. In 
the same Table we also give our values for the nucleon and A masses. As can be seen the mass of 
the A carries the largest errors. 



Volume 


(™ u ,rf) sea 


(am s ) sea 


(am q ) DW 


m™ (GeV) 


m„/m p 


m N (GeV) 


m A (GeV) 


20 3 x 32 


0.03 


0.05 


0.0478 


0.606(2) 


0.588(7) 


1.392(9) 


1.662(21) 


20 3 x 32 


0.02 


0.05 


0.0313 


0.502(4) 


0.530(11) 


1.255(19) 


1.586 (36) 


28 3 x 32 


0.01 


0.05 


0.0138 


0.364(1) 


0.387(7) 


1.196(25) 


1.643(63) 



Table 1: Mass parameters for the sea (Asqtad) and valence (DW) quarks and corresponding meson masses 
taken from [ph. In the last two columns we give our values for the nucleon and A mass using 125 configura- 
tions for am u ^ = 0.03, 75 for am u ^ = 0.02 and 38 for am u j = 0.01. 



3. Results and outlook 



quenched, K=0.1554 



f 4 



Q 2 (GeV 2 ) 



GM, 



1 



m u =0.03, m s =0.05 (125confs) 
m u =0.02, m s =0.05 ( 75 confs) 
m u =0.01, m s =0.05 ( 38 confs) 
quenched, chiral limit 



1.0 
Q 2 (GeV 2 ; 



Figure 1: G M \ in units of (e/2m N ) in Figure 2: G M \ in units of {ejlmti) using 
quenched and Nf = 2 + 1 QCD at m„ ~ 0.5 MILC configurations for three quark masses and 
GeV. quenched results at the chiral limit. 

For the 20 3 x 32 volume we analyzed 125 configurations at am u j = 0.03 and 75 at am u j = 
0.02 using both Si and S2 types of interpolating fields. This enable us to extract the dipole and 
quadrupole form factors. In addition, we analyzed 38 configurations on the 28 3 x 32 at am U i i = 
0.01 lattice using just the Si -type, providing a first indication of the magnetic dipole form factor 
behavior at a lighter pion mass. We employ the local current operator \jr(x)Yn v(x), which is not 
conserved and therefore we need to multiply by the renormalization factor Zy. This is evaluated 
non-perturbatively by evaluating the elastic nucleon electric form factor at Q 2 = 0. 

In Fig. [I] we make a comparison of the magnetic dipole form factor in quenched and dynamical 
lattice QCD at a similar pion mass of about 0.5 GeV. No unquenching effects can be observed at 
this rather heavy pion mass, given the size of the statistical errors on the dynamical QCD results. 
The results on Gm\ on the available lattices are presented in Fig. || along with the quenched results 
obtained by linearly extrapolating to the chiral limit [Q]. The accuracy of the results in full QCD 
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m u =0.03, m s =0.05 (125confs) 
m u =0.02, m s =0.05 ( 75 confs) 
quenched, chiral limit 



H 



CMR 



1.0 

Q 2 (GeV 2 ) 



m u =u.03, m s =0.05 (125 confs) 
m u =0.02, m s =0.05 ( 75 confs) 
quenched, chiral limit 




Figure 3: EMR ratio using MILC configura- Figure 4: CMR ratio using MILC configura- 
tions at am uc [ = 0.03,0.02 and quenched re- tions at am li c / = 0.03,0.02 and quenched re- 
sults at the chiral limit. suits at the chiral limit. 

do not allow for a chiral extrapolation but the general trend is that they reach lower values as 
compared to the quenched case, a welcome fact given that the quenched results were higher than 
experiment On the other hand, the EMR and CMR ratios calculated on these ensembles and 
presented in Figs. || and (| are still very noisy prohibiting a meaningful comparison to the quenched 
results. 

In conclusion, preliminary results are presented for the N to A magnetic dipole transition form 
factor Gmi and the EMR and CMR ratios, evaluated using Domain wall valence quarks and dynam- 
ical MILC configurations. The statistical fluctuations are enhanced in full QCD requiring a larger 
number of configurations in order to reach the accuracy required for a meaningful comparison to 
either the quenched data or experiment. 
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